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Abstract 

We calculate the fragmentation function for a light quark to decay into a 
lepton pair to leading order in the QCD coupling constant. In the formal 
definition of the fragmentation function, a QED phase must be included in 
the eikonal factor to guarantee QED gauge invariance. We find that the 
longitudinal polarization fraction is a decreasing function of the factorization 
scale, in accord with the intuitive expectation that the virtual photon should 
behave more and more like a real photon as the transverse momomentum of 
the fragmenting quark increases. 
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The QCD factorization theorems for inclusive single-particle production ^ guarantee 
that the differential cross section at sufficiently large transverse momentum Pt has the 
approximate scaling behavior da/dP^ ~ 1/-Pt- The scaling part of the differential cross 
section for producing the particle H with momentum P has the form 

da[H{P)\ = t dzda[i(P/z),fi]D^ H (z,fi) , (1) 

where the sum is over partons i and da is the differential cross section for producing an on- 
shell parton i with momentum Pj z. The fragmentation function gives the probability 

for a virtual parton with invariant mass of order /z to decay into a state that includes 
a particle H with the fraction z of the longitudinal momentum of the parton. The fi- 
dependence of the fragmentation functions is governed by Altarelli-Parisi evolution equations 
whose kernels can be calculated using perturbation theory. Large logarithms of z^l/Pt in 
da can be summed up by taking fi of order Pt / z. 

The factorization formula (|l[) can serve as an operational definition of the fragmentation 
functions. However, the fragmentation functions can also be given formal definitions as 
vacuum-to-vacuum matrix elements of operators that involve projections onto asymptotic 
states that include the particle H. Curci, Furmanski, and Petronzio defined them as matrix 
elements of bilocal operators in the light-cone gauge |J. Collins and Soper introduced gauge- 
invariant definitions of the fragmentation functions B. They are defined as matrix elements 
of nonlocal gauge-invariant operators that involve a path-ordered exponential of the gluon 
field called an eikonal factor: 

£ = Vexp Ug J dln^A^T^ , (2) 

where T a is a generator of the appropriate representation of the gauge group. In the light- 
cone gauge n^A^ = 0, the eikonal factor collapses to 1 and the definition reduces to that of 
Curci, Furmanski, and Petronzio. 

In the case of a light hadron H, the fragmentation function D^h cannot be calculated 
using perturbative methods. However, once Di^n(z, fio) has been measured as a function 
of z at some initial scale the Altarelli-Parisi equations can be used to evolve it to other 
scales fi. An analysis of the fragmentation functions for light hadrons has recently been used 
to obtain a high-precision determination of the QCD coupling constant Q. 

In the case of a heavy quarkonium state H, the fragmentation functions can be calculated 
using perturbation theory up to a few nonperturbative constants. The NRQCD factorization 
formalism [[J can be used to express the fragmentation function in the form 

Di-+ H (z,ii) = J2 d i^(QQ) n ( z ^)(°n) , (3) 

n 

where the sum is over color and angular momentum states of a QQ pair and the constants 
(O^) are called NRQCD matrix elements. The coefficients ^j^(QQ) n of the NRQCD matrix 
elements can be calculated using perturbation theory. The first explicit calculations were 
the coefficients of (Oi( 1 5'o)) in the fragmentation functions for g —>■ rj c and c —>■ i] c and the 
coefficients of (Oi( 3 Si)) in the fragmentation functions for g — > J/ijj and c — > J/tp ||[7|]. 
They were calculated by Braaten, Cheung, and Yuan using the operational definition of 
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the fragmentation functions provided by the factorization theorem ([[]). The first explicit 
calculations using the formal definition of the fragmentation functions were carried out 
by Ma M. The formal definition is particularly useful for calculating the fragmentation 
functions beyond leading order in a s |§. 

The factorization theorems for inclusive single-hadron production can be applied equally 
well to photon production and to lepton-pair production. They imply that at sufficiently 
large Qt, where Qt is the transverse momentum of the photon or lepton pair, the differential 
cross section has the approximate scaling behavior do jdQ T ~ 1/Q T . The scaling part of the 
differential cross section has the form ([]]). At leading order in the QED interaction, only a 
photon term needs to be added to the sum over partons i. The fragmentation functions D q ^ 
and D g ^ for real photons are nonperturbative, and only their evolution with /i is calculable 
in perturbation theory. The same is true of the fragmentation functions for a lepton pair with 
invariant mass Q comparable to the scale A associated with nonperturbative effects in QCD. 
On the other hand, if Q is much larger than A, the fragmentation functions D q ^£+£- and 
D g ^i+e- are completely calculable using perturbative QCD. Qiu and Zhang have recently 
introduced formal definitions for "virtual photon fragmentation functions" D q ^* and D g ^* 



that are equivalent to fragmentation functions for lepton pairs [T(J . 

In this paper, we calculate the fragmentation function D q ^e+£- for a quark to decay into 
a lepton pair with large invariant mass at leading order in QCD perturbation theory. We 
calculate the fragmentation function using both dimensional regularization and an upper 
limit on the invariant mass of the fragmenting quark. We find that the fragmentation func- 
tion defined by dimensional regularization has unphysical behavior except at asymptotic 
values of the factorization scale. The fragmentation function defined by the upper limit on 
the invariant mass gives a longitudinal polarization fraction that decreases as the factoriza- 
tion scale increases, in accord with the intuitive expectation that the virtual photon should 
behave more and more like a real photon as the transverse momentum of the fragmenting 
parton increases. 

The fragmentation function for q — > £ + £~ can be calculated using the operational defi- 
nition provided by the factorization theorem, as in Refs. ||[7|. The Feynman diagrams for 
producing a lepton pair with large transverse momentum Qt include some with the topol- 
ogy shown in Fig. [I], where the blob represents the parts of the diagram that involve the 
production of the virtual quark. The terms in the matrix element that corresponding to the 
fragmentation of a light quark with electric charge e q have the form 

^u{k'){ie q ej^T ^ 

where k, k', qi, q 2 , and Q = qi + q 2 are the momenta of the virtual quark, the final-state 
quark, the positive lepton, the negative lepton, and the lepton pair, respectively. The factor 
T is the 4x1 Dirac spinor associated with the blob in Fig. [j], from which the virtual 
quark emerges. If we choose a covariant gauge for QED, the numerator factor in the photon 
propagator reduces to = g^. However, if we square the matrix element (£|) and integrate 
over phase space, isolating those terms with the proper scaling behavior at large Qt, we will 
not get the correct fragmentation function. The reason is that in a covariant gauge, there 
are diagrams that do not have the topology of Fig. [TJ, but still contribute to the scaling part 
of the cross section. In order to extract the complete scaling contribution, it is necessary to 
use a light-cone gauge in which the numerator of the photon propagator has the form 



M = —u{k'){ie q ej^T—-^u{q 1 ){-iej 1/ )v(q 2 ), (4) 
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where n is the light-like vector that defines the longitudinal momentum fraction of the £ + £~ 
pair: z = n ■ Q/n ■ k. Note that the n^Q u term in (|5|) does not contribute to the matrix 
element (Q), because the leptons are on-shell. However, the Q^,n u term does contribute, 
because the decaying quark is virtual. The calculations of the fragmentation functions for 
c — > T] c , J /ip in Ref. 0] involved diagrams with the topology of Fig. [TJ, except they had a 
virtual gluon instead of a virtual photon. The authors calculated the diagram using the 
light-cone gauge propagator for the gluon. In a calculation using the formal definition, the 
contributions corresponding to the additional terms in the light-cone gauge propagator come 
from diagrams in which the virtual gluon is emitted by a gluon field from the eikonal factor 
(0)- 

If the fragmentation function for a lepton pair is calculated using the formal definition, it 
is necessary to generalize the QCD eikonal factor in (^) to include a phase from the photon 
field. In the fragmentation function for a quark with electric charge e q , the eikonal factor 
becomes 

£ = V exp (ig J dl n fl A*T a + ie q e J dl n^A^j . (6) 

At leading order in a s , the q — > C + i~ fragmentation function is given by the square of the 
sum of the Feynman diagrams in Fig. |^. The circles represent quark operators, and the 
double lines represent the eikonal factor. A convenient set of Feynman rules for calculat- 
ing fragmentation functions is given in Ref. ||. The contributions that correspond to the 
additional term in the light-cone gauge propagator (|5|) for the photon are provided by the 
diagram in Fig. [2]b in which the virtual photon is emitted by a photon field from the eikonal 
factor (H). 

In the calculation of the fragmentation function for q — * £ + £~, the integral over 
the relative momenta of the lepton pair for fixed total momentum Q gives the tensor 
~9 fiu + Q^Q V /Q 2 - In Ref. ||10|| , Qiu and Zhang decomposed this tensor into two terms 
that correspond to the polarization tensors for virtual photons that are transversely and 
longitudinally polarized with respect to the direction of the momentum of the lepton pair 
in a particular frame: 

~ 9 H Tvj — [~9 H ^ ) + ™ ^ ^ {') 



Q 2 V a ' 2 J \ Q 



where n is a conjugate light-like vector satisfying n ■ n = 2. The frame chosen in Ref. |10 
is the one in which the perpendicular momenta of the decaying quark and the lepton pair 
are k' ± = —Q±/z and Q' ± = 0. In this frame, the second tensor on the right side of (0) can 
be expressed as e^e£, where the polarization vector is proportional to (n ■ Q) 2 n^ — Q 2 n^. 
This choice for the decomposition into transverse and longitudinal contributions is rather 



artificial, but we will adopt it in order to compare our results with those of Ref. [JTO . 

We proceed to calculate the fragmentation functions from the square of the sum of 
the two Feynman diagrams in Fig. |2|. We use dimensional regularization in 4 — 2e space- 
time dimensions to regularize ultraviolet divergences. The fragmentation function can be 
expressed as an integral over the invariant mass Q 2 of the lepton pair: 
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,<0) , Ce]a 2 i dQ 2 (4n^ 2 \ 2 ' \ 2(1 - z) + (1 - e) 



z 2 



^«-w = -sN-y # (# j { ' v w.w - ami 



2 2 ) 

where C is a constant that reduces to 1 as e — > 0: 

(l-2 £ )(l-f £ )r(i-e)' ^ 

In (§), there is an implicit lower limit on Q 2 that is large enough that a perturbative 
calculation of the QCD corrections would be reliable. The two terms in the integrand 
correspond to transversely and longitudinally polarized photons, respectively. The functions 
f n {z) can be expressed as integrals over the perpendicular components Qx of the momentum 
of the lepton pair. Equivalently, they can be expressed as integrals over the invariant mass 
s of the decaying quark, which is related to by light-cone energy conservation: 

z l-z V ; 

This implies the lower limit s > Q 2 /z. The functions f n {z) in (§) are given by 

"(1 - e) Ji/z y n 



where y = s/Q 2 . They can be evaluated analytically: 

r(w-l + e) 



fJz) = -± , r ' Z n -\l -Z)' € . (12) 

The function f\(z) in @ has a pole in e. This ultraviolet divergence is cancelled by the 
renormalization of the composite operator in the definition of the fragmentation function. 
The renormalized fragmentation function in the MS scheme is 

e l a r 1 <h 

2n Jz y 



D^ e+i .(z^) = D®(z) - f 1 d -l Pq ^ z /y)D^ l+l -(y) , (13) 



where P q ^ is the splitting function 

P„ = I±<Iz£>I, (M) 
and D^£+£- is the fragmentation function for a photon to decay into a lepton pair: 



3tt 7 Q 2 \ Q 2 



The coefficient C is given in (||). Our final result for the MS 1 fragmentation function is 
obtained by taking the limit e — » in ([13]) : 
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[Z,fi) 



6tt 2 



Q 2 



i + (i 
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z)Q 2 





+ 2(1-^)1 











(16) 



The two terms in the integrand correspond to transversely and longitudinally polarized 
photons, respectively. 

In Ref. ||10|| ) Qiu and Zhang defined their fragmentation functions by imposing the con- 
straint s < jJ? F on the invariant mass of the decaying quark, which sets the upper limit 
y < fip/Q 2 on the integral in fllTl) . This eliminates ultraviolet divergences, so we can set 
e = 0. The results for the integrals are then 




- 1 6{& - Q 2 /z) 



2 904 - Q l /z) 



- Q 2 /z) . 

The resulting expression for the fragmentation function is 

dQ 2 nr. .2 ^2,., fl + (l-^) 2 



(17) 
(18) 
(19) 



D 



e 2 g a 2 
Gtt 2 



Q 2 



e0& 



Q 2 /z) 




(20) 



This result has been confirmed by Qiu and Zhang [10]. If we take the formal limit fi 2 F 3> Q 2 / z 
in (PD]), the 9 function becomes 1 and the resulting expression differs from flTED only in the 
argument of the logarithm and in the coefficient of the z term. The arguments of the loga- 
rithm are the same if we make the identification fj/p = fi 2 / (z(l — z)). This can be understood 
by examining the expression (|T0"| ) for the invariant mass of the decaying virtual quark. It 
suggests that the scale \x of dimensional regularization should be identified not with a cutoff 
on the invariant mass of the virtual quark, but with a cutoff on the perpendicular momentum 
Q± of the lepton pair. The difference between the coefficients of the z term comes from the 
ezf\(z) term in (|8J), which reduces to z in the limit e — » 0. On the other hand, if we impose 
an invariant- mass cutoff fip and set e = 0, this term vanishes. Taking the limit fip ^ Q with 
z fixed in corresponds simply to an alternative ultraviolet cutoff. The fragmentation 
function should then differ from ( |T6"D by a finite renormalization of the composite operator 
in the definition of the fragmentation function. This finite renormalization corresponds to 
adding the term — ez to the splitting function P q -,-y in (pf). 

One advantage of the fragmentation function fl20[) defined by an upper limit on the 
invariant mass is that it builds in threshold effects and the constraints of energy conservation 
associated with the decay of a virtual quark with invariant mass This might be useful for 
quantitative applications of the fragmentation function. In Fig. |3|, we compare the differential 
fragmentation functions Q 2 dD(z)/dQ 2 corresponding to ( p0|) with invariant mass cutoff fip 



and ( |T6"D with renormalization scale \v 



1 — z)fj,p. We set e q = +^ and a 



We 
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-1 ^3 clliu u 137' 

choose Q = 5 GeV and consider 2 values of [ip: ftp — 10 GeV in Fig. |a and /ip = 50 
GeV in Fig. |3|b. The fragmentation function defined by s < fip is below the threshold 
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at z = Q 2 1 fip and is positive for larger values of z. The fragmentation function defined 
by dimensional regularization has unphysical negative values over part of its range as can 
be seen in Fig. |3|. As z — > 1, it remains positive only if fi F > exp(l)Q. As z — > 0, it 
is negative for all //p, diverging like \ia(zfip/Q 2 )/z. For large enough values of hf/Q as in 
Fig. [3|b, the two fragmentation functions look similar except that one vanishes for z < Q 2 / fi 2 F 
and the other becomes negative in that region. If we keep \i fixed in (|TJ|) instead of fip, the 
fragmentation function still exhibits unphysical behavior. It is positive-definite if \x > Q, but 
it diverges like ln(/i 2 / Q 2 ) / z as z — > and like ln(/i 2 /((l — z)Q 2 )) as z — > 1. We conclude that 
the fragmentation function ( [TB|) defined by dimensional regularization is of little practical 
use. It is essential to take into account threshold effects in some way, such as by imposing 
an upper limit on the invariant mass as in (P0[). 

If the QED phase in the eikonal factor (|6]) is omitted in the formal definition of the 
fragmentation function, D^ e+e _ is given by the square of the Feynman diagram in Fig. 0a. 
The resulting expression is not gauge invariant, but it is independent of the gauge parameter 
for covariant gauges. Using dimensional regularization, we obtain 



*^-^> ~ 6tt 2 J Q 2 \ Q 2 



+ 



z 

l-Z 2 



2z 5 



z 2 f (z)-Azh(z)+Af 2 (z) . (21) 



The two terms in the integrand correspond to transversely and longitudinally polarized 
photons, respectively. The transverse term is identical to that in (|16|) . This follows from 
the fact that the Feynman rule for the emission of a virtual photon from the eikonal line in 
Fig. |||b is proportional to n M , which is orthogonal to the transverse tensor in (0). Thus the 
omission of the diagram in Fig. ^b can only affect the longitudinal term. With dimensional 
regularization, the function fo(z) vanishes, as is evident from fll2~|) . The function fi(z) has 
a pole in e. In the transverse term, the pole can be removed by the renormalization (|13|). 
However, there is also a pole in the longitudinal term that is not removed by renormalization 
of the composite operator. This failure of renormalization is the signal that the definition of 
the fragmentation function that omits the QED phase in the eikonal factor is inconsistent. 

In Fig. f|, we compare the transverse and longitudinal contributions to the fragmentation 
function fl2"0"|) calculated using an invariant-mass cutoff. We choose Q = 5 GeV and consider 
2 values of fip: \ip = 10 GeV in Fig. f|a and \ip = 50 GeV in Fig. £|b. The dashed curves 
labelled T and L are the transverse and longitudinal contributions given by the two terms 
in (P0|). Their sum is the solid curve. The longitudinal polarization dominates just above 
the threshold at z = Q 2 //! 2 ?, because the longitudinal term increases linearly in z — Q 2 /fi 2 F 
while the transverse term increases quadratically. The transverse polarization dominates at 
large z, because the longitudinal term vanishes at z = 1. 

We can define a longitudinal polarization fraction £ by dividing the longitudinal fragmen- 
tation probability by the total fragmentation probability for lepton pairs of invariant mass 
Q 2 . The fragmentation probability is proportional to the integral over z of dD(z)/dQ 2 . The 
longitudinal polarization fraction for Q = 5 GeV decreases from £ = 0.56 at fip = 10 GeV 
to 0.34 at fi F = 50 GeV and to 0.23 at fi F = 250 GeV. This is in accord with the intuition 
that the virtual photon should behave more and more like a real photon as fxp increases. 

The longitudinal polarization fraction defined by the ratio of fragmentation probabilities 



7 



decreases rather slowly with fi F . However a more relevant measure of the polarization is 
the ratio of the second moments of the fragmentation functions. The reason is that the 
hard-scattering cross section da scales like 1/(Qt/z) 2 , and this weights the fragmentation 
function by z 2 . The longitudinal polarization fraction for Q = 5 GeV defined by the ratio 
of the second moments decreases from £ = 0.43 at fip = 10 GeV to 0.15 at fip = 50 GeV 
and to 0.08 at y. F = 250 GeV. 

An important observable in lepton pair production is the angular distribution of the 
momentum of one of the leptons in the rest frame of the lepton pair. The angular distribution 
is proportional to 1+a cos 2 9, where 9 is the angle between the momentum qx of the negative 
lepton and some quantization axis. The polarization variable a is related to the fraction £ of 
lepton pairs that come from longitudinally polarized virtual photons by a = (1 — 3£)/(l + £). 
The fraction £ depends on the choice of quantization axis. The choice adopted in Ref. |T(J is 
rather artificial, because it requires specifying the transverse momentum of the fragmenting 
quark which is not easily observed. In the case of a hadron collider, a more physical choice 
for the quantization axis is the direction of the momentum Q of the lepton pair in the rest 
frame of the colliding hadrons. In this case, the longitudinal polarization vector is a linear 
combination of Q and the total momentum K of the colliding hadrons: 

Q 2 K»-(K-Q)Q» 



[(#.Q) 2 -i^Q2]l/2[Q2]l/2- 



The contribution to the fragmentation function from longitudinally polarized virtual photons 
can be obtained by replacing the lepton tensor — g^ v + Q fl Q u /Q 2 by e^e u L . The resulting 
expression for the fragmentation function depends explicitly on k + = k ■ n: 

e 2 a 2 r dQ 2 2(1 — z) 

D q ^ (e+i - )L (z,fi F ) = -j-^- J — — f 2 (z,fx F /Q,Q/k + ) , (23) 

where ji is a function of z, fip/Q, and Q/k + : 

f U a IO Oik ) - r 1 ' * dy " VMML (24 ) 



In the limit k + 3> Q, this reduces to (|T9|). The complete fragmentation function summed 
over polarizations is independent of Q/k + and is given by (p0[): 



In the frame defined by K = 0, k + is the sum of the energy and momentum of an on- 
shell quark produced by some hard scattering. A reasonable choice for \x F is the transverse 
momentum of that quark, give or take a factor of 2. 

The NRQCD factorization approach predicts that the 1 quarkonium states should be- 
come increasingly transversely polarized as their transverse momentum P F increases HTT| , |T2"11 • 
Quantitative predictions of the polarization of the ip{2S) [0,0, J/ip jl5l| , and T(2S) [16| 
indicate that the increase in the polarization should set in at values of P F that are accessible 
at the Tevatron. The present data on the polarization of J/if> and ip{2S) from the CDF 
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collaboration |0| seem to indicate a decrease in the transverse polarization at large Pt, 



although in both cases the discrepancy with the prediction is significant only in the largest 
Pt bin. The argument that the tranverse polarization of J/ip or ip{2S) should increase with 
Pt is completely analogous to the corresponding argument for lepton pairs, except that 
it involves a virtual gluon instead of a virtual photon. There are many effects that could 
dilute the transverse polarization or delay the onset of the predicted increase. However no 
plausible mechanisms have been identified that could make it decrease with Pt- We expect 
that more accurate measurements from Run II of the Tevatron will reveal the increase in 
transverse polarization predicted by NRQCD. 

In conclusion, we have calculated the fragmentation function for a light quark to decay 
into a lepton pair to leading order in a s . For renormalizability and for QED gauge invariance, 
it is essential to include a QED phase in the eikonal factor in the formal definition of the 



fragmentation function. Berger, Gordon, and Klasen [18[ have shown that the distribution 
of the transverse momentum Qt of lepton pairs in hadron collisions is dominated by parton 
processes initiated by gluons if Qt > Q/2. The Qt distribution can therefore provide useful 
constraints on the parton distribution for gluons. Our fragmentation function for q — > £ + £~ 
may be useful for calculating the Qt distribution in the limit Qt S> Q- 
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FIG. 1. The class of Feynman diagrams that in the light-cone gauge for QED gives the 
contribution to the production of lepton pairs associated with fragmentation of the light quark q 
at leading order in a s . 
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FIG. 3. The differential fragmentation function Q 2 dD(z) / dQ 2 for q — > for a lepton pair 
with invariant mass Q = 5 GeV as a function of z for (a) fj,p = 10 GeV and (b) fj,p = 50 GeV. 
The solid curves are the fragmentation functions ([H]) defined by the invariant-mass cutoff /J,p. 
The dashed curves are the fragmentation functions ( |i~6| ) defined by dimensional regularization with 
renormalization scale /x 2 = z(l — z)[i 2 F . 
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(a) (b) 

FIG. 4. The differential fragmentation function Q 2 dD(z)/dQ 2 for q — > for a lepton pair 
with invariant mass Q = 5 GeV as a function of 2 for the invariant-mass cutoffs (a) fip = 10 GeV 
and (b) = 50 GeV. The solid curve is the total fragmentation function. The dashed curves 
labelled T and L are the contributions from transverse and longitudinal virtual photons. 
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